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Abstract 

The SL(2,R) WZNW — * Liouville reduction leads to a nontrivial phase 
space on the classical level both in + 1 and 1 + 1 dimensions. To study 
the consequences in the quantum theory, the quantum mechanics of the + 1 
dimensional, point particle version of the constrained WZNW model is in- 
vestigated. The spectrum and the eigenfunctions of the obtained — rather 
nontrivial — theory are given, and the physical connection between the pieces 
of the reduced configuration space is discussed in all the possible cases of the 
constraint parameters. 



1 Introduction 

In the past several years the Toda models have been studied intensively. In these 
field theories scalar fields are coupled to each other by certain special exponential 
terms, in a way that corresponds to a simple Lie algebra. The Toda models can 
be considered as generalizations of the Liouville theory, which is of particular inter- 
est since it appears in many problems of physics and mathematics. An interesting 
means of deriving and studying the remarkable properties — integrability, conformal 
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invariance and W-algebraic symmetry — of the Toda models is offered by the ob- 
servation that the Toda theories can be obtained by a suitable reduction of the 
Wess-Zumino-Novikov-Witten (WZNW) model Q. The WZNW model is a theory 
of a field that takes its values from a Lie group G, and the reduction procedure - by 
imposing appropriate (first class) constraints - associates it to a Toda theory that 
corresponds to the Lie algebra of G. In the case G = SL(2, R) the reduced theory 
is nothing but the Liouville theory. 

However, the connection between the WZNW and Toda models is more intri- 
cate. A closer look at the reduction procedure shows that it yields not exactly the 
Toda theory but a richer structure, the Toda theory arises only as a component, a 
subsystem of it. This aspect was first noticed in [|]. To study the precise relation of 
the Toda models to the WZNW ones a recent work examined the SL(2, R) WZNW 
— > Liouville reduction from the phase space point of view ||. The authors con- 
sidered the classical SL(2, R) WZNW model, imposed the appropriate constraints 
and described the reduction of the phase space under the constraints. They found 
that the reduced phase space contains two subsystems (nonintersecting open sub- 
manifolds) that admit a clear physical interpretation. On both subsets the reduced 
WZNW theory leads to the Liouville theory locally, but these two copies of Liouville 
theories are not independent. The connection between them comes from a 'border 
line', a lower dimensional surface in the reduced phase space connecting them. For 
a better understanding of the situation || carried out a similar analysis on the + 1 
dimensional, point mechanical analogue of the SL(2, R) WZNW model. This can 
be thought as the space independent version, the 'zero mode sector' of the WZNW 
model. In [[3] it was found that depending on the signs of the constraint parameters 
one can arrive at two different types of reduced theories. In both cases the phase 
space reduces into two locally independent parts. The difference is that when the 
constraint parameters have equal signs the two halves are disconnected, there is no 
'border line' between them. As a result a classical motion cannot touch both parts. 
Actually, though the reduced Hamiltonian is not of the usual form of the sum of a 
kinetic and a potential term, the system behaves as if the two halves of the reduced 
one dimensional configuration space were separated by an infinitely high potential 
barrier. On the other hand, if the signs of the constraint parameters are opposite 
then an analogue of an infinitely deep potential well characterises the situation. In 
this case the two 'half worlds' are connected, a motion can cross the point that 
separates the two halves of the reduced configuration space, moreover, the negative 
energy motions will oscillate between the two parts. Q concludes that it is not 
enough to give the (global) reduced theory in the local coordinates of the reduced 
phase space. The arising two components of the reduced theory seem to be inde- 
pendent while actually they may have a physical connection, a correlated behaviour 
which can be discovered only from the global point of view. (|| also contains results 
about the SL(n, R) point particle model and discovers some similar properties for 
n > 2 as well.) 

These are the features of the classical theories. What happens on the quan- 



2 



turn level? Because of the expected difficulties of quantizing the SL(2, R) WZNW 
model it is useful to examine the point particle version. Turning to the masspoint 
theory the following questions arise naturally: what kind of relation connects the 
two half-worlds in the quantum theory? Are there any oscillating motions quantum 
mechanically when the the constraint parameters are of opposite sign? Similarly we 
may ask if in the case of coinciding signs the separation remains or a tunneling is 
allowed. Moreover, in the first case one may expect negative energy bound states. 
Are they really present? In this paper we solve the quantum theory of the masspoint 
version, and thus we can answer whether our expectations based on the classical be- 
haviour hold. While — as a consequence of its rather nontrivial characteristics — the 
point particle problem is interesting in itself, its properties are expected to shed a 
light on the field theoretical case, just as it happened on the classical level ||. 

The paper is organized as follows. Sect. 2 and 3 present the classical mechanics 
of the SL(2, R) masspoint before resp. after imposing the point mechanical form of 
the WZNW — > Liouville constraints. The quantum mechanics of the unconstrained 
and the constrained systems are established in Sect. 4 and 5. The reduced system 
splits into two parts in a symmetric way, Sect. 6 gives the eigenf unctions on one 
such part. Sect. 7 investigates the orthogonality and the completeness properties of 
these 'half-eigenfunctions'. The connection between the two parts is examined in 
Sect. 8. Sect. 9 discusses the coordinate independence of the definition of the con- 
strained system. Sect. 10 gives the conclusions of the paper. The larger proofs and 
calculations belonging to the statements of Sect. 7 are presented in two appendices, 
one concerning orthogonality and one concerning completeness. 



2 The classical mechanics of the unconstrained 
theory 

The SL(2, R) WZNW model is defined by the following action: 



where a and r coordinate a two dimensional Minkowski space, g is an SL(2, R)- 
valued function of o and r being periodic in a with period 2ir. The coefficients of 
the first term, the action of the SL(2, R) sigma model, and the second, topological 
term called Wess-Zumino term, are denoted by m/An and m! respectively. The point 
particle version, i.e. the physics of the zero modes of the field theory ([!]) is defined 
by restricting the configurations to the space- independent ones g = g(r) only. Then 
the action reduces to 



Tfl f r i 

5o = — ) dadri:i[{g- l d T gf - (g'^g) 2 ] + 






(2) 
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The motion of the point particle is a function g : R — > SL(2,R); we can see that 
the Wess-Zumino term does not contribute to the masspoint version. The theory 
posesses left and right translation symmetries under the transformations g ^ hg, 
g i— > ghr x , h e SL(2,R), the corresponding conserved quantities are J = gg^ 1 and 
J := g~ x g, taking their values in si (2, R), the Lie algebra of SL(2, R). 

The equation of motion following from the action (0) is (g~ l g)' = (or, equiva- 
lently, {gg~ x )' = 0). Its solution is g(t) = g(0)e At , where A e sl(2,R) is a kind of 
initial data specified by the initial conditions as A = g(0) _1 g(0). Using the basis 
T\ = o\ y T 2 = i(T2, T 3 = o"3 in sl(2, R) (where <7fc-s denote the Pauli-matrices) A can 
be written as A k Tk (Einstein-summation). We will make use of the identity 

e cfe ^- = cosh R1 + sinh R/R c k a k (3) 

where R 2 = |Tr[(c fc <7fc) 2 ] (this formula also holds for Tr[(c fc <Tfc) 2 ] < with R being 
imaginary and for Tr[(c fc <Tfc) 2 ] = with coshi? — > l,sinhi?/_R — > 1). Based on @ 
the solution of the equation of motion is 

/ ,\ _ / n \ ( cosh (rt) + A 3 sinh (rt)/r (A 1 + A 2 ) sinh(rt)/r \ , . 

9\t) ~ g(u) y ( A i _ A 2^ sinh ( rt ) j r cos h(rt) - A 3 sinh(rt) jr ) ^ > 

with r 2 = |Trp fc T fc ) 2 ] (also allowing zero or imaginary values of r). If the trace is 
negative then the solution also can be expressed as 

n(A - „(n\ ( cos( ^ + A3 sin M/^ ( Al + A 2 ) M"t)/u \ 

9W ~ 9{V) I ( A i _ A 2^ gin ^ cos ^ _ A 3 gin ^ t ) I W 



with lo = \J —r 2 . For r = the solution actually reads 

U\ 1 + AH (A 1 + A 2 )t \ . 

g(t) = g(o) [ {A i _ A 2 )t X _ AH J- (6) 

A compact form of formula (f|) is g(t) = g(0) (cosh(rt) 1 + sinh(r£)/r A). We can 
see that the r 2 < solutions are closed because of their trigonometrical time depen- 
dence, the r 2 > motions are open. 

To study the canonical structure of the theory let us consider a parametrization 
of SL(2, R) g(C)i i = 1, 2, 3. The Lagrangian reads as 

£({,0 = ^(0^, (7) 

where hki{£) = Tr [g~ 1 dkgg~ 1 dig] is the metric tensor on SL(2,R). The canonical 
momenta are 

dL 

Pk = T7T = mh kl£ l , (8) 

from which we get the Hamiltonian 

H(^p) = ^-h kl p kPl (9) 
2m 
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({h kl }, the inverse of the matrix {hki} exists because SL(2, R) is simple). 

The Hamiltonian can be expressed as the function of £ and £ as well. One finds 
that 

m 

Combining this result with the coordinate free form of the Lagrangian (cf. eq. (fj)) 
and g{t) = g(0) e At , the value of the Hamiltonian on a solution of the equation of 
motion can be determined easily: 

H = yTr[v4 2 ] = mr 2 . (11) 

As a consequence we can see that the energy is negative for a closed motion and 
nonnegative for open motions. 

Let us turn to a special parametrization, namely the one which is based on the 
Gauss decomposition: 



a (3 \ ( 1 a \ ( 1/5 \( 1 
7<fJ~loiHo 5 I [ c 1 



(12) 



This parametrization describes any y ^ ^ J e SL(2, R) except those having 5 = 0. 
With 5 = £\ a = £ 2 and c = £ 3 {/ifcj is of the form 

{h k i} 

The determinant of this metric is 



/ 2/5 2 

5 2 | . (13) 
V 5 2 



/i = -25 2 , (14) 

and the inverse of the metric tensor, {h kl } has the form 

/ 5 2 /2 \ 
{h kl } = 1/5 2 . (15) 
V 1/5 2 / 

Thus we arrive at the following Hamiltonian in this parametrization: 

n(5,a,c,ps,p a ,Pc) = j-5 2 p 2 + (16) 

Expression (|i~6| ) shows that a and c are cyclic coordinates since 7i is independent of 
them. This is an advantage of using the parameters 5, a, c. Later we will see that 
this parametrization fits very well for our further considerations. That's why in the 
following we will work in these coordinates. 
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3 The constrained model on the classical level 



Now we impose the constraints that reduce the SL(2, R) WZNW theory to the 
Liouville one. For space independent configurations they read 



Tr 



engg 



Tr 



6213 l g 



;i7) 



where eia : 
tion gives 




and e 2 i 




. Expressing (|17D in the 5, a, c parametriza- 



p a = mfi, 



Pc 



mv. 



As p a and p c are constants of the unrestricted motion we see that the constraints 
mean nothing else but a special choice of some of the initial conditions. This is a 
general feature of first class constraints, that these constants of motion happen to 
be canonical momenta is the real advantage of the parameters 5, a, c. 

The reduced phase space can be obtained by factorizing the complete phase space 
by the gauge transformations these first class constraints generate. The action of 
these gauge transformations is 
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(19) 



(see ||) where 9i and 9r are the two parameters of the transformations. This 
symmetry transformation acts on a and c as a i— > a + 9l, c h c + 9r, while it 
leaves 5, ps, p a and p c invariant. Thus the factorization simply means that S and ps 
parametrize the reduced phase space. 

Any motion S(t),ps(t) allowed by the constrained dynamics corresponding to the 
constraint parameters n and v can be obtained by taking a solution 8(t),a(t),c(t), 
Ps(t),p a ,p c of the unconstrained theory where p a = m/i, p c = mv. Roughly speaking 
we just have to omit a{t) and c(t). The coordinate 5 G R survives the reduction, thus 
one may consider it as the coordinate of the one dimensional configuration space of 
the reduced theory. The constrained dynamics is governed by the Hamiltonian 



(20) 



where s = \iv. 

Until now we have found the Gauss decomposition a very appropriate way of 
introducing coordinates on the group SL(2, R) to reach the canonical structure of 
the reduced system. Now let us face the problematic side of this approach. 

The topology of SL(2,R) is not R 3 but R 2 x S 1 , consequently it cannot be 
covered by a single parametrization. In the case of the Gauss decomposition the 
signal of this is that the Gauss decomposition works only for 5^0, the elements 
( a (3 s 
[ 7 



are left out. As a result the canonical coordinates 6, a, c, p$, p a and p c 
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parametrize only two nonintersecting open submanifolds, two 'open halves' of the 
whole phase space, corresponding to the two regions — oo < 5 < and < 5 < oo. 
In particular, the point 5 = is left out from the reduced configuration space. In 
the meantime, extracted from (|j), the time dependence of 5 is of the form 5(t) = 
C\ sinh(rt)/r+C2 cosh(rt) where C\ and C*2 are arbitrary constants. Thus we can see 
that there exist motions that cross the 5 = surface in the whole phase space — for 
example for imaginary values of r oscillations occur between the regions 5 > and 
5 < 0. Furthermore, for such a motion ps tends to infinity as the particle reaches 
5 = 0. This can be seen from ps(t) = 2m5/5 2 , the connection between p$ and the 
smoothly varying quantity 5 (cf. (§) and (|13|)). Hence the coordinates 5, p$ of the 
reduced phase space seem to be able to describe only those parts of a motion when 
5(t) ^ 0, they cannot give an account of how the particle moves across 5 = 0. 

Fortunately we can overcome these difficulties. First, let us complete the reduced 

configuration space by mapping the points ^ a ^ ^ of the configuration space of 

the unconstrained system to the point 5 = of the reduced configuration space, 

just as we have mapped the points ( ^ ^ ^ ^ ^ to 5. Second, let the recipe to 

follow a motion through 5 = be to fit the quantity 5 2 ps- This recipe is clear 
from the unconstrained point of view: we simply fit 5 here. By these tricks we 
can eliminate the disadvantage of working with only one parametrization instead of 
covering SL(2, R) with more than one patches. Clearly, the canonical formalism of 
the unconstrained system also owes the problem of the missing 5 = 0. Nevertheless, 
the solution given here for the reduced system applies in a straightforward way for 
the unconstrained one, too. 

We will see that the 5 = problem also arises in the quantum theory. There it 
appears as an irregular singularity of the Hamiltonian at 5 = and the challenge 
will be to define the quantum theory on the whole configuration space despite this 
singularity. 

Finally let us introduce a canonical transformation which transforms the con- 
strained Hamiltonian to a form of a sum of a kinetic and a potential term. This can 
be achieved by the following transformation: 



x 



:=V2\n5, (21) 



p x := -±=5p 5 . (22) 



The resulting Hamiltonian is 



H x (x,p x ) = ^p\ + mse-^ x . (23) 

The price we have to pay for having such a nice Hamiltonian is that by (^T|) we 
restricted ourselves to 5 > only. (Or, because of the 5 «-> —5 symmetry of the 
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system, to 5 < 0, if writing —5 instead of 5 in fl2ip.) Remarkably, the logarithmic 
connection between x and 5 is the point mechanical analogue of the one that relates 
the field of the reduced SL(2,R) WZNW theory to the Liouville field <fi in the field 
theoretical case 0. 

With the aid of fl23|) it is easy to analyse the three qualitatively different situa- 
tions arising. If s > then the potential increases exponentially as we travel to the 
negative x direction. Thus for all the allowed motions with positive energies there 
is a turning point when moving to the negative direction towards x = —oo (5 = 0). 
In this case there is no possibility for the masspoint to cross the border 5 = 0. For 
s = we have a free particle moving along the x axis. Now the 'point' x = — oo 
(5 = 0) cannot be reached in a finite time interval so the masspoint cannot cross the 
border even in this case. In the case s < an exponentially deep potential valley 
attracts the particle towards the negative direction, what's more, the time needed 
to reach x = — oo happens to be finite. This shows that for s < the particle may 
cross the border 5 = 0. 



4 The quantum theory of the unconstrained sys- 
tem 

Let us define the quantum mechanics of the point particle SL(2, R) WZNW theory 
via canonical quantization. We use the coordinates 5, a, c and work in the coordinate 
representation. The wave functions are then complex valued functions defined for 
all 5 G R \ {0}, a G R and c G R. We define the scalar product as 

(V u ifr 2 ) : = J ^ 2 V^hd5dadc. (24) 

Here h denotes the determinant of the matrix {hki} in the 5, a, c parametrization 
(cf. (|14"D). The measure in this integral is the usual one used on curved manifolds. 
Moreover, in our case the metric tensor is invariant under the left and right transfor- 
mations. Consequently, if we adopt the left and right transformations for the wave 
functions: 

[D L (h)*](g) := ^(h^g), [D R (h)*](g) := *{gh), g,he SL(2, R), (25) 

then the scalar product is also invariant. This property is inevitable if we want 
the left and right symmetries of the classical theory to be present on the quantum 
level as well. Observe that these natural requirements led to the appearance of a 
nontrivial weight function p(5) := \J — h = V2\5\ in the integral (p4l). 

In coordinate representation the canonical momenta are defined as partial deriva- 
tions with respect to the corresponding coordinates if the configuration space is flat. 
For curved configuration spaces this definition does not give symmetric operators 
because of the presence of the weight function in the scalar product. Symmetricity 
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and the requirement [ik,Pi] = ifrbki meet in the definition pk : = j{dk + \dk In \f—K) 
(see || for example). In our case this gives p§ '■— ^(dg + ^r), p a := hd a and p c := ~d c . 

Now the stage is set for the definition of the Hamiltonian. This is also a step 
which needs some care. The problem with the naive quantum analogue of the 
classical Hamiltonian ([16]): 

Un = ^~S 2 p 2 s + -PaPci- (26) 

is that 7i n is not symmetric (do not forget the weight function!). We need to find an 
appropriate ordering of 5 and p$. Moreover, symmetricity is not the only requirement 
since we would like the Hamiltonian to be invariant under the left-right symmetry 
transformations as well. That is why we turn to the Laplacian of SL(2, R): 

A *=vb|j(^ ( «^*)' (27) 

the usual definition for curved manifolds, which reads in the case of the 5, a, c 
parametrization 

A ^ = W\ dsm3 m) + i dadc ^- (28) 



The definition ( 27|) and the properties of the metric tensor guarantee that the Lapla- 



cian is invariant not only under reparametrizations but under left and right trans- 
formations as well. The reason why the Laplacian is interesting is that — |j^A is 
just a re-ordered version of 7i n . So we define the Hamiltonian as 

-Ml) (5,a,c) := -^-l- d5 f\ S \^^ aiC )) _ ^ _L d a d c V(5, a, c). (29) 
' Am | o | v J m o z 

A simple calculation shows that this Ti is symmetric (with respect to the weight 
function p). This choice of 7i is the usual and natural one for curved configuration 
spaces (see, for example, [|J). 

We mention that we have not defined the wave functions at 5 = and the 
operators 5, ps and H are apparently ill-defined at 5 = 0. These singularities are 
only coordinate artifacts here. This won't be the case for the reduced system as we 
will see soon. 



5 The quantum mechanics of the reduced system 

Being ready with the quantum mechanics of the unconstrained theory the next task 
is to consider the quantum analogue of the constraints and see what the reduction 
yields. Let us impose the constraints on the quantum level as 

Pa%! = m[i^, p c ^ = mv^ (30) 
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(cf. flTB|)). It is very easy to find the wave functions that satisfy (0), they are of 
the form 

*(<J, a, c) = ip(5) e ^ a+muc \ (31) 

Remarkably, the action of the operators S, p$ and TC on a wave function of this form 
touches only its 5-depending part. This makes it possible to work with i/j(5) instead 
of ty(5,a,c) and to consider the one dimensional quantum mechanics driven by a 
Hamiltonian H: 

m)(5) = -^^Hl^dsH^+rns^S) (32) 
together with the scalar product 

(Vi,^):= [r^2P(S)d5 (33) 



Therefore in the following we investigate the properties of this one dimensional 
quantum system. (The classical H corresponding to this H is just ( p0|) as we expect.) 
Here we witness how the decoupling of the variables a, c from the system happens 
on the quantum level. 

The wave functions ( |3~1~1) are not square integrable in the SL(2, R) sense. This 
is a natural consequence of the constraints that decrease the degrees of freedom by 
two, we will require square integrability 'in the reduced sense', i.e. with respect to 
the scalar product (|33|). The situation is similar to the case of a free masspoint 
in a three dimensional Euclidean space with the constraints classically imposed as 
p y = const., p z = const. When we consider the corresponding quantum theory it is 
obvious that the normalizability of the wave functions must be understood 'in the 
one dimensional sense'. 

The one dimensional problem we arrived at is quite an unusual one. The Hamil- 
tonian is not a 

ij> ^ -b 2 ip" + Vip (34) 

-type (with a real constant b and a potential function V) and a nontrivial weight 
function is present in the scalar product. It would be very convenient if our system 



could be transformed to an 'ordinary' one with a Hamiltonian of the form (|34|) and 
with no weight function. For this purpose let us consider a transformation of 5 to a 
new variable x: 

x = g(5) (35) 

accompanied by a change of ip(8) to a new wave function xi x ) which is related to 

as 

m = f(S)x(g(5)), (36) 

where / and g are arbitrary real functions. To see how H transforms under such 
a transformation one has to examine the transformation of the (time dependent) 
Schrddinger equation : 

zh^(6,t) = (H^)(6,t). (37) 
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After a simple calculation the result can be written in the form 

ift^OM) = -^(grx"-^-(s6g' + 26^g' + 6 2 g")x' + 



dt 4m 4m y / ' 

X- (3* 



4m 1 / / / <5 2 



Requiring eq. (p8f) to have the form given by eq. (p4|) gives a number of equations. 
First of all the coefficient of the second derivative must be a negative constant, for 
convenience we want to set it to — j-- This yields 

S 2 (g') 2 = 2. (39) 

The solution of this condition is 

g(5) = V2\n( Cl 5) (40) 

where C\ is an arbitrary positive constant; for the moment let us restrict ourselves 
to the positive 5 half line. Secondly we require that the x' term must vanish. This 
and fl40|) together yield 

f(S) = j (41) 



with a nonzero constant 0%. Under the transformation fl36|) with such an f(S) and 
g(S) the weight function changes as 

52 ft(6)i/> 2 (6)V2\6\d6 = f 9(52) xUx) X2 (x)cldx. (42) 

We want to transform away the weight function entirely. This can be reached simply 
by choosing c 2 = 1. The parameter c\ does not have such an interesting specific 
value, we set C\ = \. The transformation of the variable 



x 



\/21n5 (43) 



and the wave function 



lead to the following H x : 



m = \x(x) = \x(V2\n5) (44) 



(H xX )(x) = -^X"(x) +(^ + mse-^A x (x). (45) 



This transformation is just the quantum analogue of (^T|) and (23). (The ad- 
ditional constant in the potential term of H x is due to the ordering procedure we 
maintained at the definition of the quantum Hamiltonian.) Unfortunately the prob- 
lem is the same as well: it works only for the positive half of the configuration space 
(or the negative one, if exploiting the symmetry 5 <-> —5). Nevertheless, H x will be 
very useful to understand the physics encoded in H. 
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6 Eigenfunctions on the half configuration space 



With the aid of the transformed Hamiltonian ( f45[) one can have a rough picture of 
the reduced theory. In the cases s > 0, s < it is more or less similar to a system 
with a potential infinitely increasing or decreasing for 5^0, while for s = the 
system is somehow a 'sum' of two free theories. The potential valley of the case 
s < suggests to use the Bohr-Sommerfeld quantization to get a first impression 
about the spectrum of H. For this purpose let's consider a classical motion with 
energy E and express ps as a function of S using (f20|): 

V '4mE5 2 — 4m 2 s , , 

ps = ± j 2 . (46) 

For the Bohr-Sommerfeld quantization condition the phase space area § p$d5 is 
needed as a function of E. The problem is that for 5 ~ the behaviour of ps is 
Ps ~ 1/5 2 (here and in the following « means asymptotical or approximate equality 
and ~ means proportionality), so the integral is infinite. Thus the Bohr-Sommerfeld 
quantization is impossible. 

What is the reason behind this? The answer is that H is not bounded from 
below if s < 0. To see this let us consider a square integrable wave function if) and 
define 

iP x (S) := \iP(\S). (47) 
The if)\-s are normalized to 1; the expectation value of H in a state if)\ is 

(V>A, Hlpx) = (VA, + C0Aj ^2^A), (48) 

where f/i and ff 2 denote the first and second term of H, respectively (cf. (|32|)). 
The scaling properties of Hi and H 2 are such that 

(ifj x ,H^ x ) = ^,H^) (49) 

and 

/oo 1 
\iP(5)\ 2 -V2\S\d8. (50) 
-co 

In the r.h.s. of (|50D A 2 is multiplied by a negative number. As a result if A increases 
to oo then (ip\,Hip\) tends to — oo. 

This lack of a ground state causes the failure of the Bohr-Sommerfeld quan- 
tization. Fortunately the eigenvalues and eigenfunctions of H can be determined 
exactly in all the cases s > 0, s = 0, s < 0. To do this we have to solve the equa- 
tion Htf) = Eif) as a differential equation of second order. This equation has three 
singular points: 5 = ±oo, which are regular singular points and 5 = 0, which is an 
irregular sigular point. Consequently one has to solve this equation restricted to 
the domains 5 G R + and 5 G R~ respectively and then to fit together the obtained 
'half-eigenfunctions'. Because of the symmetry 5 <-> —5 it is enough to work on R + , 
the restriction of H to R + will be denoted by H + . 
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In the case s < the eigenvalue equation can be transformed to the Bessel 
equation 

z 2 w" + zw' + (z 2 - v 2 )w = (51) 



by the substitutions 



and 



where 



and 



z = ~ § (52) 



m = ]*(£\, (53) 



v 2 = 1 - (54) 



^f^. (55) 

The two linearly independent solutions of fl5l|) — existing for any complex value of 

z/ — are J u (z) and 3^(z) (for the conventions and properties concerning the Bessel 
functions cf. ||). Similarly, for s > the transformations 

z = g, (56) 

m = \v(j) (57) 
lead to the modified Bessel equation 

z 2 w" + zw' - (z 2 + v 2 )w = (58) 

with v 2 being the same as in (|54]) and 



/ 4m 2 s 

*=y-^-- (59) 

Now the two solutions of ([58]) are the modified Bessel functions I u {z) and In 
the case s = the transformation is the most useful. The eigenf unctions of H x 



are exp(±iKx), where K = \jlmE j% — 1/2 (cf. (|45|)). In the variable 5 they read 

1 (60) 

In the cases s > 0, s = if^ is bounded from below. The corresponding 
condition on the energy eigenvalues is E > j-. For s = it means that K is a 
nonnegative real number. In the case s > the condition gives v 2 < causing that 
only the functions Ii U (z) and K iu (z), u E R mean energy eigenf unctions. For s < 
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the Bessel functions with real indexes lead to E < j- and the imaginary indexes 
correspond to E > 4-. 

To get more acquainted with the eigenf unctions let us carry out a simple check 
of our physical picture that is based on H x . In the cases s > 0, s < the potential 
term of H x decreases exponentially to zero as x tends to oo. Consequently we expect 
that for x — > oo the eigenfunctions with E > j- behave as plane waves. (For s = 
this expectation is satisfied trivially.) To see whether this is the case we will make 
use of the z ~ behaviour of Ji u (z) and Ii U {z). Using ( |96|) up to a constant of 
proportionality they are of the form 

J iu {z) » I lu {z) ~ (61) 

From fl43|) and ( |52|) the connection between the variables x and z is 

z = fce"^, x = V21b^, (62) 

(or, for s > 0, the similar formulas with k instead of k,) so we can see that in the 
variable x Ji U and Ii U are asymptotically plane waves. The momentum correspond- 
ing to them is p = —hu/\^2. Considering that for x — > oo V(x) tends not to zero 
but to Ti 2 /Am and quoting the connection between v = iu and E we find that the 
expectation 'kinetic energy = p 2 /2m' is satisfied as well. The two other eigenfunc- 
tions, Yi U and Ki U are linear combinations of Ji U and J-i U resp. 1{ U and I- lu . Thus 
they also behave the way we expect from our physical picture. 



7 Orthogonality and completeness 

It will be important to form a complete orthogonal system from the half-eigenfunctions, 
an orthogonal basis in L 2 (R + , p). In the case s = this is simple: the set {exp(±iKx) \ K e 
[0, oo)} is a complete orthogonal system (in the variable x 6 (—00,0x3)). Therefore 
the same can be said about the functions (|60|) in the variable 5 in L 2 (R + ,p). For 
s < it is shown in the appendices that there exist several independent choices of 
a complete orthogonal system. The different bases can be indexed by a p e (0,2], 
the corresponding eigenvectors (given in the transformed form w(z)) are 

J q (z), q = p, p + 2, p + 4, 

cos (fp) J Q (z) + sin (fp) Y (z), (63) 
e^ w 4(«) + ^ w /-»(4 u E (0, 00), 



where 



el e p {u) = cos (f g) sinh(f m) + i sin(f p) cosh(f u) 

/cos 2 (|p) sinh 2 (|-u) + sin 2 (|p) cosh 2 (|w) 
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For the case s > the appendices prove that only one complete orthogonal system 
can be built from the functions Ii U (z) and K iu (z), namely, the set 



{K iu {z) | u G [0,oo)} 



(65) 



What makes the difference that in the cases s = and s > the eigenbasis is 
unique while for s < there are infinitely many complete orthogonal systems? The 
answer is in the self-adjointness of H + . For this reason we determine the deficiency 
index of H + . H + is a differential operator of second order with real coefficients and 
two singular points 5 = 0, 5 = 00. Its deficiency index is equal to the number of its 
orthogonal square integrable eigenfunctions corresponding to a non-real eigenvalue 
(cf. ||) (the deficiency index does not depend on the eigenvalue chosen). In the 
case s > the deficiency index is zero as for a fixed non-real v 2 none of the two 
linearly independent eigenfunctions — I v {z) and K v (z) in the variable z — is square 
integrable (cf. Appendix A). In the case s < J u is square integrable while Y v is not 
(we can choose Rei^ > without loss of generality, see Appendix A again). Thus in 
this case the deficiency index is one. For s = the deficiency index is zero, which 
can be seen most easily in the variable x. 

Now, starting with the case s < 0, we recall a theorem of 0, which states that 
if the deficiency index is one then the operator has several self-adjoint extensions. 
also gives a condition for the different domains of definition of the different self- 
adjoint extensions. For H + this condition says that a function ip{5) lying in the 
domain of definition of a self-adjoint extension has (to be smooth enough, cf. ||, 
and) to satisfy 



lim 

<5^0 




d8 



^ 7/ 



lim 



~dJ 



where 



(66) 



(67) 



with a $ G [0, 2tt) and a v G C \ R, Re v > 0. $ and v together index the different 
self-adjoint extensions. 

Then if one examines which eigenfunctions are included in the domain of defi- 
nition of a self-adjoint extension indexed by an arbitrarily chosen value of 1? and v, 
a straightforward if lengthy calculation shows that these eigenfunctions are exactly 
the ones that form one of the complete orthogonal systems flBlf). The number p 
which characterises this system is expressed by d and v as 



p = Re v 



arcsm 



sinh(|Im v) sin 



7T 



'sinh 2 (f Im v) sin 2 (f ) + cosh 2 (f Im v) cos 2 (|) 



(mod 2) 
(68) 

(For deriving this result one can make use of the asymptotics of the Bessel functions 
and their derivatives, for the asymptotics cf. Appendix A). Hence the multiplicity 
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of the eigenbases origins in the multiple self-adjoint extensions of the differential 
operator H + . 

In the cases s > and s = the deficiency index is zero. The appropriate 
theorem of states that then the operator is self-adjoint. Consequently, the domain 
of definition is unique. All the eigenfunctions (or, more precisely, all the wave packets 
superposed from the eigenfunctions — remember that for both s > and s = all 
the eigenfunctions are non-normalizable) are lying in the domain of definition so the 
eigenbasis is unique as well (up to linear equivalence, i.e. except from trivial phase 
factors or, in the case s = 0, choosing two linear combinations of exp(iKx) and 
exp(— iKx) instead of exp(iKx) and exp(— iKx)). 

8 The eigenfunctions on the whole configuration 
space 

To investigate the eigenfunctions of the full system the task is to sew together the 
half-eigenfunctions and build up a complete orthogonal system of 'whole-eigenfunctions' 
(in the following: eigenfunctions). In usual quantum mechanical systems, i.e. with 
a Hamiltonian of the form ( p4|) and with no weight function in the scalar product, 
the conditions for fitting parts of an eigenfunction together are the continuity of the 
eigenf unction and the continuity — or in special cases a given jump — of its (space) 
derivative. Now we cannot expect that such conditions work. In fact, the 5 — > 
behaviour of the half-eigenfunctions proves to be cos(k/5 + const.) in the case 
s < 0, 5 _1 exp(±i\/2K In 5) if s = and 5~z exp(— k/8) if s > 0. Thus this kind of 
fitting together is impossible. The situation is not better in the variable x either, 
the half-eigenfunctions tend to in the limit x — > — oo in the cases s < and s > 0, 
while for s = they behave as exp(±zi^a;). This infinite growth or decrease and in- 
finitely rapid oscillating behaviour of the half-eigenfunctions origins in the irregular 
singularity of the Hamiltonian at § = 0. 

Fortunately, the probability current is finite at 5 — >• 0, it is this quantity we are 
able to fit. However, in our case the probability current is not of the usual form 
though. By deriving the continuity equation for the probability density from the 
Schrodinger equation the probability current proves to be 



Any ip can be expressed as a linear combination of eigenfunctions (p^, which makes 
it possible to decompose the probability current as a sum of 



It can be verified that in each cases s<0, s = 0, s>0 such a quantity has a well 
defined finite limit for 5 —>■ so the probability current is also finite at 5 — > 0. 




(69) 




(70) 
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We do not fit the probability current directly but carry out an equivalent proce- 
dure. In fact, fitting the probability current of the half-eigenfunctions is to ensure 
that the norm of a whole wave function does not change in time. The latter is 
equvivalent to the self-adjointness of the whole Hamiltonian. We know that the 
eigenf unctions of a self-adjoint Hamiltonian are orthogonal. Conversely, a complete 
orthogonal system of the eigenfunctions of the Hamiltonian as a differential operator 
defines a self-adjoint Hamiltonian from the differential operator on an everywhere 
dense set in L 2 (R,p), which is our purpose. That's why it is enough — while it is 
more interesting as well — to build up complete orthogonal systems out of the eigen- 
functions instead of fitting the probability current. 

Let us start with the case s < 0. An eigenfunction <&(5) is generally of the form 

J a<p(-6) if 5<0, 

\ pip{5) if (5 > 0, 1 1 

where ip is a half-eigenfunction, defined on R + . From this it follows immediately that 
at most two linearly independent eigenfunctions can correspond to an eigenvalue in 
a (whole-)eigenbasis. Another important observation is that if a value p corresponds 
to (p — the index of the half-eigenbasis (p is a member of — , this p characterizes $ as 
well. Now let us suppose that a complete orthogonal system of eigenfunctions does 
not include two linearly independent eigenfunctions that correspond to the same 
eigenvalue and have the same value p. (Later we will examine the other well, 
i.e. when one can find two such eigenfunctions in the system.) In this case there 
must be at least one eigenfunction in this eigenbasis with a different p. Otherwise 
we do not have completeness: there exist functions that are orthogonal to any basis 
vector but are not identically zero; such an example is a whole-eigenfunction that is 
not included in the basis but has the same value p. 

For two eigenfunctions having different p-s ( , )+ and ( , )_, the restriction of 
their scalar product to the positive resp. negative half of the configuration space 
are not zero, consequently they are orthogonal only if one of them is of the form 

const { if 5 < 0, 

COnst -\ \<p(5) if 5>0, [U) 

and the other is of the form 

f -X*ip(-5) if 5 < 0, /f7nS 

(with a different ip but) with the same complex A from the set {|A| < 1, if |A| = 
1 then argA G [0, 7r)}. Let p\ denote the value p of the eigenfunction of the first form 
and P2 the p of the other one. The other eigenfunctions with p\ also must have the 
form ( |72|) (with the same A) and the other eigenfunctions with p<i also must have the 
form (173]) , in order to be orthogonal to these two eigenfunctions. These forms ensure 
that the further eigenfunctions are orthogonal to each other as well. Orthogonality 
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also excludes the existence of any eigenfunctions in the eigenbasis having a p other 
than pi or p 2 . 

To examine completeness first let us see whether an arbitrary function ipi from 
L 2 (R,p) having the form fl72|) (where now (p is not a half-eigenfunction but an 
arbitrary half-function) can be spanned by these eigenfunctions. It is easy to see that 
this requirement is equivalent to that the restriction of the eigenfunctions with p\ to 
R~ have to form a complete half-eigenbasis (a ipi is orthogonal to the eigenfunctions 
with p 2 , hence only the eigenfunctions with p\ contribute to it). After a similar 
treatment of the ip 2 s of the form (^) we conclude that a complete system must 
consist of each of the eigenfunctions with p\ (that have the form (|72|)) and each of 
the eigenfunctions with p 2 . Then if any ip e L 2 (R, p) can be given as a sum of a ipi 
and a ij) 2 then completeness is reached. With the notation 



the sum of the functions 



and 



*>o (75) 

* (i):= l ^(-a<M*)+^M), *>o (76) 

is ip, thus the completeness of the considered system of eigenfunctions — which we 
shall denote by (pi,P2,ty — is proven. We remark that the above decomposition of 
ijj is a generalization of the decomposition of a function to a sum of an even and an 
odd function, which is actually the case A = 1. 

Now let us turn to the other case, i.e. when the complete orthogonal system of 
whole-eigenfunctions includes two linearly independent eigenfunctions with a same 
value p and corresponding to the same eigenvalue. In this case the other eigenfunc- 
tions must be of this p as well, otherwise they cannot be orthogonal to both of these 
eigenfunctions. Furthermore, the eigenvalues of the eigenbasis must be identical 
with the eigenvalues of the half-eigenbasis p and must be doubly degenerated: in 
the case of a simply degenerated or missing eigenvalue any (other) eigenfunction cor- 
responding to this eigenvalue is orthogonal to each eigenfunction from the system, 
which is in contradiction with completeness. The constants a and (3 (see (fFTD ) for 
the eigenfunctions of the system can be arbitrary, the only requirement is that for 
each eigenvalue the corresponding two eigenfunctions be linearly independent. The 
concrete values of these as and j3s are not important, they only embody a choice 
of two basis vectors in a two dimensional linear subspace. Remarkably, such an 
eigenbasis is linearly equivalent to a one which — in the spirit of our notation — can 
be denoted by (p,p, A) (the equivalence holds for an arbitrary A). Based on this 
observation one can prove completeness the same way as for a system (pi,p 2 , A). 
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We see that in contrast to the eigenbases (pi,P2, A) considered earlier, these latter 
eigenbases are characterized by a single number p. Each of the different eigenbases 
(pi,P2,X) and p means a different self-adjoint extention of the Hamiltonian as a 
differential operator. 

For s > and s = the method to establish an eigenbasis is the same as for an 
eigenbasis with a single p in the case s < 0. The difference is that now one starts 
with only one half-eigenbasis, consequently one arrives at only one eigenbasis (up to 
linear equivalence). As a result in these cases the self-adjoint Hamiltonian is unique. 

That for s < two different half-eigenbases are needed in general for one eigen- 
basis may seem unusual. However, this situation is just an analogue of the case of 
the operator —d 2 /dx 2 on the interval [— 7r, it]. This simple example will help us to 
understand - at least to some extent - what happens here. 

The complete orthogonal system of the eigenvectors of —d 2 /dx 2 corresponding 
to the conditions 

$(-tt) = $(tt) = (77) 



consists of the functions 



sin 



n . . 

-(X + Vr) 



n = l,2,... (78) 



These eigenfunctions are either even or odd functions, the even ones can be written 

as 

sinf-z], fc = 2,4, (79) 



2 



the odd ones are of the form 




/ = 1,3,... (80) 

Let us examine how one can get this eigenbasis by building it up from the half- 
eigenbases of —d 2 /dx 2 , the eigenbases of the operator restricted to [0, n] and [— %, 0]. 
The eigenfunctions of [0, tt] corresponding to the conditions 

<p(Q) = <p(ir) = (81) 



arc 



m = l,2,... (82) 



We see that these functions can be the building blocks of the functions (|79|), via an 
antisymmetric (i.e. odd) extension from [0, it] to [— ir, tt]. On the other side, there 
is no way to build up the functions ( jSOD from them as well. The solution is that 
we have to consider another half-eigenbasis, namely, the one corresponding to the 
conditions 

y/(0) = <p(n) = 0. (83) 
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The eigenfunctions satisfying (|83| ) are 

1 3 

cosjx, j = -,-,..., (84) 
and clearly a symmetric (i.e. even) extension of them to [— 7r, it] yields the functions 

The operators — d 2 / dx 2 and H are common in that both commute with the space 
reflection operator, and in that both have several self-adjoint extensions. These are 
the aspects why the operator —d 2 /dx 2 is a good example to understand H. However, 
crucial differences can be observed between the two operators. One of them is that 
the proof that a Hamiltonian with space reflection symmetry must have even and 
odd eigenfunctions does not hold for H. The other, more serious difference is that in 
the case of —d 2 /dx 2 there exists a physical principle to choose one from the different 
self-adjoint extensions. This is possible because the operator —d 2 /dx 2 together with 
the interval [— n, n] arises in physics as the Hamiltonian of the system characterized 
by the potential 

V{x) = {° (85) 

v I OO if \X\ > 7T. 

Then the requirements that an eigenfunction be continous and that it be zero where 
V(x) = oo assign one of the possible boundary conditions that characterize the 
different self-adjoint extensions (this distinguished condition is just (ff7|)). In this 
case there is a way to choose the 'physical' self-adjoint extension of the Hamiltonian. 
For H there is no such principle, all the self-adjoint extensions prove to be equal. 
We have to face the fact that there is no unique quantum mechanics corresponding 
to the s < classical system. 

In spite of the non-usual form of the Hamiltonian and the presence of the nontriv- 
ial weight function, the H x form of the Hamiltonian enables us to give the physical 
interpretation of the results to some extent. 

In the case s < we expect that the two half-configuration spaces are in physical 
connection, the particle can cross the border 5 = 0. Simple calculations show that 
this expectation is satisfied for the self-adjoint extensions (pi,p 2 , A), there is a prob- 
ability flow from one half to the other one. Consider for example a wave function 
ip = C\ipi + C21P2 where ip\ and ip2 are eigenfunctions, one having pi and the other 
having p 2 . Though f t {ip,ip) = 0, f t (ifj,ip) + = -f t (ip,ip)- ^ for generic ci and 
C2. Another transparent possibility to show the physical connectedness of the two 
halves is that one can easily find examples for a solution of the (time dependent) 
Schrodinger equation where the expectation value of the coordinate operator 5 is 
oscillating in time between a positive and a negative value. However, in the case of 
the self-adjoint extensions p the two halves behave as two closed, independent sub- 
systems. The reason is that for these eigenbases the restriction of the eigenfunctions 
on a half configuration space is a half-eigenbasis, causing that H decouples to two 
self-adjoint half-operators. 
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The cases s > and s = are similar to the s < 0, p one. The Hamiltonian is 
simply a pair of two self-adjoint half-Hamiltonians, the two parts of the configuration 
space are physically independent. This result is in accord with the naive pictures 
of the s > and s = systems based on H x . For s > we can think of an 
exponentially increasing and thus infinitely wide potential wall separating the two 
half-worlds, no wonder that we find no tunneling from one side to the other. The 
situation is similar to the quantum mechanics of the system with the potential 

V(x) = [° ^h >a ' (86) 
v ' [ oo if \x\ < a. v ' 

where it is meaningful to speak about the quantum mechanics of the system on 
the whole configuration space, yet there is no physical connection between the two 
allowed parts. In the case s = we have two free theories, both on an infinitely 
large configuration space (understood in the variable x). We may argue that under 
such circumstances a wave packet starting from one side (e.g. the 5 > one) cannot 
reach the other side in a finite time period. We cannot say any stronger concerning 
interpretation: these are the limits we are forced into. 



9 Coordinate independence 

The reduced quantum system, as we saw, posesses several unusual properties. It 
is natural to ask whether these features are only artifacts, caused by the special 
coordinate system which was used for the definition of the reduced system. Therefore 
it is worth examining the possibility to define the system in a coordinate independent 
way. 

The Hamiltonian of the unconstrained quantum theory — a multiple of the Lapla- 
cian of the manifold SL(2,R) — and the scalar product ( [24|) are in fact coordinate 
invariant. Consequently the question reduces to whether the constraints can be given 
a coordinate independent form. In Sect. 5 the constraints were imposed through the 
canonical momentum operators. The definition pk '■— j (dk + hdj. In y/—h) does not 
define a covariant quantity because h is not a coordinate invariant scalar. That's 
why it is recommended to impose the constraints independently of the canonical 
momentum operators. 

In the spirit of the Lie derivative, let us introduce the following derivation oper- 
ators 



for any A G sl(2,R). The definition of La and Ra does not need any coordinate 
system. Nevertheless, if expressing them using the coordinates 5, a, c one finds that 
L ei2 = d a and R e21 = d c . Thus we obtained a coordinate independent reformulation 
of the constraints (BO). 
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One can feel the need for checking whether the operators j^L ei2 and j^R e21 are 
really the quantum equivalents of the classical quantities Tr [en gg~ x ] and Tr [e2i g^g] 
The following heuristic argument makes this relation visible. 

We consider a wave packet which is in some sense the most similar to a classical 
trajectory, determine the expectation value of ■A-L^, and compare it to Trf/l^g -1 ] 
computed on the classical trajectory the wave packet is similar to (the relation will 
hold for any A £ sl(2, R) in general; the analogous treatment for Ra is straightfor- 
ward) . 

In the case of the free quantum mechanics on a three dimensional Euclidean 
space the Gaussian wave packet 

const. J ^ke'^'^e^-^ (88) 

is in some sense the best wave mechanical analogue of a classical trajectory. It is 
well-localized both in position and in momentum — around the position x and the 
wave vector k$. We wish to define the analogue of this wave function in the case 
of SL{2, R). The differences in the two configuration spaces are that SL(2, R) does 
not have a linear structure and the natural metric on it is not positive definite. 
Nevertheless, we 'flatten' SL(2, R) around a chosen point go — which we want as the 
'centre' of the wave packet — by characterizing a point g with a Lie algebra element 
denoted by X g defined by the relation g = e Xg go. This definition is correct, i.e. 
X g exists and is unique, if g is in a small enough neighbourhood of go, and this 
is enough for our purposes because we want to define only a fairly well localized 
wave packet around g . X g is the analogue of x — x in the Euclidean case. As 
far as the indefinite metric is concerned: that a wave packet is localized is actually 
a question involving topology rather than metric. That's why we will not use the 
natural metric {hki}, or, more precisely, the corresponding Killing form G, to express 
that our wave packet is localized. (G is the metric tensor at the unity element of 
SL(2,R), acting on the Lie algebra of SL(2,R), G(U,V) = Tr[UV] ). Instead, 
we introduce an arbitrary positive definite inner product G in si (2, R) and apply 
this G in the Gaussian modulus of the integrand in (|88|). This procedure is similar 
to the case of the Minkowski space, where for example differentiability or other 
analytic properties are defined not by the use of the Lorentzian metric but by an 
arbitrarily chosen positive definite inner product; it turns out that these properties 
are independent of the choice of this artificial inner product. In the phase of the 
integrand of (|88|) we keep the Killing form. 
All in all, let us work with the wave packet 

V go>Ko (g) = const. J d 3 Ke-^ K - K °> K - K ^ e iG(K ^l (89) 

The width a will not have to be specified. This wave packet is by definition localized 
around K , evaluating the integral shows that the result is a Gaussian in position 
(also in the sense of G), localized around go as expected. 
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The function g i— > exp iG(K, X g ) is approximately an eigenfunction of -^La with 
eigenvalue — G(K, A). By this we mean that the eigenvalue equation holds exactly 
at g = go and is satisfied approximately at a neighbourhood of go. Consequently the 
expectation value of j^La is —G(Ko, A) — approximately, i.e. in the limit when the 
width of the wave packet in position tends to zero. Then what is left is to show that 
—Kq = gg^ 1 for the classical trajectory corresponding to (|89|). To obtain this we 
make use of the observation that A = (G~ 1 ) l - ? Lj i Li j = (G _1 ) lJ R^Ri- where is 
an arbitrary basis in sl(2, R). The most straightforward verification of this formula 
is to prove it in coordinates. From this one can derive that g i— > expiG(K, X g ) is 
approximately an eigenfunction of H = —^A with eigenvalue ^G(K, K). With 
the aid of these pieces of information we are able to tell the time propagation of the 
wave function (|89|) . If considering (j39|) as the state function at t — then after a 
small time period the state function is approximately 

* 90jKo (g,t) = const. J#Ke-&* K - Ko *- K °)<W K *h-i$**™). (90) 

is also a localized wave packet, but its centre is not go but a g(t). This g(t) 
plays the role of the classical position of the masspoint at t. Using a stationary 
phase argument and taking into account that the modulus of the integrand of ( |90"D 
is concentrated around K the centre of the wave packet ( ]9ll| ) is at X g ^ = ^K t. 
Consequently, gg~ l = —K , and this is what we wanted to show. 

We saw that both the unconstrained system and the constraints are actually 
coordinate independent. After defining a system in a coordinate independent way, 
one can use concrete parametrizations to examine its properties. Turning to the 
concrete situation: all the properties explored in the coordinates 5, a, c are valid 
everywhere where 5^0. For example, the wave functions are scalars so the infinite 
growth and infinitely rapid oscillating of the eigenfunctions is a coordinate indepen- 
dent fact, since this is the behaviour of the eigenfunctions not at but around the 
invalid point 5 = 0. 

10 Conclusions 

We investigated the properties of the point particle version of the reduced SL(2, R) 
WZNW model both on the classical and the quantum level, for all the possible values 
of the constraint parameters. We found that the quantum theory exhibits an analo- 
gous behaviour to the classical one. The cases where the two parts are disconnected 
classically lead to two independent systems on the quantum level as well, and in the 
cases where the half-systems have a physical connection, this connection can also 
be found in the quantum theory. The only exception is that there is a possibility 
for a classically connected case to be disconnected quantum mechanically. This is 
possible because not only one quantum theory corresponds to a classically connected 
case. Several self-adjoint extensions of the Hamiltonian exist, including special ones 
where the two half-systems turn out to be independent. 
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Classical mechanically bounded motions exist, with arbitrary large negative en- 
ergies, in the connected cases. The disconnected cases do not allow bounded mo- 
tions and energy is bounded from below. These properties are also reflected on 
the quantum level. It is remarkable that the quantum theory is formally consistent 
irrespective of the values of the constraint parameters, while in the connected cases 
it leads to systems with a Hamiltonian not bounded from below (no matter which 
self-adjoint extension is chosen). Recently a method was proposed to discuss quan- 
tum mechanical systems that exhibit such a behaviour [0] . The method implements 
the concept of Wilson renormalization. It would be interesting to carry out such an 
analysis for the system studied here. Nevertheless, the method of J7| means a kind 
of distortion of the system, which is not the purpose here as here we are interested 
in the properties of the original system for we want to obtain indications how the 
quantum theory of the corresponding field theory behaves. In Sect. 3 we have found 
classical space-independent configurations with arbitrary large negative energy in 
the connected cases. This and the quantum properties of the masspoint version 
make it quite possible that the energy is essentially not bounded from below in the 
quantum field theory. 

That the energy is not bounded from below is not the only nontrivial property 
of the connected case. The most striking result of our analysis is the existence of 
several self-adjoint extensions corresponding to one classical system. There is no - 
physical or mathematical - principle to choose one out of them as the 'real' one. 
The origin of this behaviour is the strong singularity at the border which separates 
the two half-systems. This singularity is not present on the unconstrained level, it is 
a consequence of the charasteristics of the constraints. As this singularity can also 
be observed in the classical reduced field theory version B, we expect to face the 
problem of the non-unique self-adjoint energy operator on the quantum level, i.e. in 
the quantum field theory of the reduced SL(2, R) WZNW model as well. 

The method applied here to present the quantum mechanics of the reduced 
system was canonical quantization (supplemented by a coordinate independent ap- 
proach). Because of the nontrivial properties found it would be interesting to exam- 
ine this system by using other tools, geometric quantization or functional integration, 
and see how these methods give account of the characteristics of the theory. 

Additionally we remark that recently a paper carried out an analysis of the 
relativistic quantum mechanics of a free particle on the SL(2,R) manifold ||. The 
problem studied there is independent from the one presented here. Clearly, in 
the group SL(2,R) plays the role of the (curved) space-time the particle exists in 
while in our case SL(2, R) is the (configuration) space of the unconstrained system. 
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A Orthogonality 

Let us first collect the properties of the Bessel functions and modified Bessel func- 
tions we will need in the following. J u (z) and Y u (z) are analytic functions of z on the 
whole complex plane except the negative real half line. For a fixed z (z ^ 0) both 
are integer functions of v for any v G C. J u (z) and J- U (z) are linearly independent 
for any value of v except v = n (let n denote an integer subsequently), when 

J_ n (z) = (-1)" J n {z). (91) 

J v {z) and Y v (z) are always linearly independent. The following formula 

= ™(™)Uz)-J-u(z) 

establishes a connection between the J-s and the Y-s; for v = n ( |9~2"D is understood 
as a limit when v — > n. 

For the modified Bessel functions I v {z) and K u (z) all the statements mentioned 
above apply if we replace J u (z) by I v (z) and Y v (z) by K u (z), with the exceptions 
that I- n (z) = I n (z) and that 

M2) = !M!hW. (93) 

2 sm{nu) 

Under complex conjugation 

J u (z)* = J l/ *(z) (94) 

is valid if z is real. This property holds for Y v (z), I v {z) and K v {z) as well. 

The modified Bessel functions can be expressed in terms of the ordinary Bessel 
functions. We will make use of the identity 

I v (z)=e~ i ^J u {e i ^z) (-7r<argz< |). (95) 

The asymptotic behaviour of these four functions at z ~ is 



J v {z) « I v {z) « — — y , (i/ G C \ {-1, -2, -3, ...}), (96) 
Y v (z)n--K u (z)a-~r{v)(~) \ (Rez/>0) (97) 



IT 7T \2 

and 



F (^) « --K (z) w -In 2. (9? 

7T 7T 
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For \z\ — > oo the ordinary and modified Bessel functions behave as 



JJz) 



Y„'(z) ps \ — cos 



Jv'(z) 



-Y v {z) « 
h(z) 



TTZ 

2 



7T 




sm 



7TZ 



71 



(|argz| < tt) 
(|argz| < tt), 



v/2 



7T 

arg^l < - 



and 



~ -Kj{z) 



are; 2 < 



3tt 



(99) 

(100) 
(101) 

(102) 



Expression ( P5[ ) is valid not only for fixed z/ and 2; ps but holds also if z is fixed 
at an arbitrary (not necessarily small) value and \v\ — > 00, |argz/| < tt. In this case 
it can be combined with the Stirling formula giving the \u\ — > 00 asymptotics of the 
gamma function resulting 



J„(z) PS I v {z) 



V2 



7T1S 



ez 
2^ 



(103) 



After these necessary pieces of information let us start finding the possible or- 
thogonal systems of the eigenvectors of the case s < 0. The scalar product of ipi(5) 
and i>2{&) in L 2 (R + ,p) is 



(^1,^2 



(104) 



here the notation ( , ) + reminds us that this scalar product is taken in L 2 (R + ,p), 
i.e. on the positive half of the configuration space only. Instead of the variable 5 
it will be more suitable to work in z. Under the transformations (p2| ) and fl5H| ) the 
integral ( 104|) transforms to 

V2 [°° w 1 (z)*w 2 (z) — . (105) 
Jo z 

We will study this integral by considering it between finite a and b and then take 
the limit a — > 0, b — > 00. 



If W\(z) and W2(z) are eigenfunctions of the Bessel equation (51]) with indexes ^ 
and respectively then such an integral can be easily evaluated due to a formula of 
HI. This formula states that 



AJz)BJz 



dz 



v 2 — fJ? 



(A,(z)Bl(z) - A^z)B v {z) 



(106) 



(A, B = J or Y, < a < b < 00). Applying it together with ( |9"41 ) the integral ( |105| ) 
between the limits a and b is 



— ^ \ z (wi(z)*w' 2 (z) - w' 1 (z)*w 2 (z))] b a . 



— [i 



(107) 
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With the aid of ( |107| ) let us determine the scalar product (J^, J u ) + (/x, v G C). 
Substituting them into (|107|) and using the asymptotic formulas (|99|) and (|100|) , the 
contribution from the 'upper limit terms' (the terms depending on b) in the limit 
b — > oo is 



2^ 



7TI V 



*2' 



sin 



7T 



:io8i 



On the other hand, from 



one can see that the 'lower limit terms' behave as a 



7/ 



-th power of a when a — > 0. If Re (//* + z/) > then the 'lower limit terms' 



tend to zero; so in this case 



2^2 



7r(z/ 2 — yU* 2 



Sill 



7T 



(109) 



If Re (/!* + 1/) < then the integral diverges for a 
limit gives the norm of J^. 



0. If Re n > then the v — > /i 



v^Re// 



(110) 



Similarly we find that ( J^, 1^)+ is finite if Re (//* + u) > and infinite if Re (p*+v) < 
0, while the scalar product (Y^, Y v ) + always diverges. As a special case of this latter 
we also see that the Y^-s are non-normalizable, for any complex value of //. 

We want to build complete orthogonal systems out of the J M -s and Y^-s, where 
now only real and imaginary /x-s are allowed (cf. Sect. 6). The role of the Bessel 
functions with real respectively imaginary index is to span the wave functions that 
have an energy expectation value E > j— resp. E < j-. With a real /j only 
the J M >o-s and Y can be taken into account. In fact, a J M <o or an Y^q cannot 
be orthogonal to a Bessel function with imaginary index as their scalar product 
diverges. From ( |109| ) it follows that a maximal, pairwise orthogonal set of J M >o-s is 
{ J p , Jp+2, J P +A, ■ ■ ■}, where p G (0, 2]. With a p given, the only linear combination 
of a Jiu(z) and J-i U {z) that is orthogonal to J p (z) is the one given in (|63D, as one 



finds with the aid of ( p,09 ). This linear combination proves to be orthogonal to the 
eigenf unctions J p +2, </ P +4, ... as well. In the end, there exists one linear combination 
of Jq(z) and Y (z) that is orthogonal to J p (and, as turns out, to the functions J p+ 2, 
J p+ 4, . . . and exp[— i6 p {u)} J iu (z) +exp[i9 p (u)} J_ iu (z), u G (0, 00) as well). This linear 
combination can be obtained either by using the scalar product or as the u — ► 
limit of the eigenfunctions exp[—i0 p (u)]J iu (z) + exp[i9 p (u)} J_ iu (z). 

What is left to check is the mutual orthogonality of the eigenfunctions with 
imaginary index. We will prove that 

(c(u) [e"*Wj m + e^J-iu] , c(v) [e~ w ^ J iv + e^W J_ fo ]) = 6(u - v) (111) 



where 



c{u) 



\T(l + iu) 



u 



2^ sinh(7TM) 



(112) 
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(for the properties of the gamma function see for example 0). As m and v run 
over the positive real numbers only and not on the whole real line, the Dirac delta 
distribution must be understood here to act on the test functions which are defined 
on the positive half of the real line and are smooth functions of compact support 
vanishing at the origin. In the following u will be treated as a variable — i.e. the 
variable of the test functions and the kernel functions — and v as a fixed parameter. 
To inspect the scalar product let us consider (|107|) in our case. A simple calcu- 



lation involving the use of ( |108| ) and some trigonometrical identities shows that the 
contribution of the upper limit tends to when b — > oo. The lower limit terms are 
asymptotically 



V2c(u)c(v) 



-i(Op(u)-0 P (v)) / a \i(u-v) 

+ ... 



(u-v)r(l + iu)T(l-iv)\2 



(113) 



where the . . . stands for three other terms which can be obtained from the first one 
by the substitutions u — * — u, v — > — v, or both, respectively. Introducing A = In ~, 
two of these terms have the form f(u, v) sin A(u + v ) + g(u, v) cos A(u + v ), where 
/ and g are smooth functions of both u and v. It is a well-known fact that the 
regular distributions sin Ax and cos Ax tend to zero if A — ► oo. This property does 
not change if we multiply them by a smooth function, thus these terms give zero in 
the limit A — > oo (a — ► 0). 

With (|11 2Q the two other terms can be written in the following way 



2m{u — v) 



/ e -i[0 p (u)-6» p («)+argr(i+iu)-argr(i+w)] _ 



-2 



-(«<-> v). (114) 



From the properties of 8 P and the gamma function it follows that the function 
[exp(. . .) — 1]/ (u— v) behaves smoothly even if u — >• v. Consequently the distributions 
coming from the first and the third terms also tend to zero. What remained is equal 
simply to sin[A(u — v)] /tt(u — v). It is well-known that sinAa;/7rx — ► S(x) in the 
limit A — > oo, so ( |1 1 1| ) is proven. 

We close Appendix A by showing that in the case s > {K iu (z) | u G [0, oo)} 
is the only possible orthogonal system built from the eigenfunctions. For this let us 
consider the scalar product 

(c + (u)I iu + c^(u)I_ iu ,c + (v)I iv + c^(v)I_ iv ) + (115) 

(here u, v > again). In the case s < we had formula ( |106| ) to evaluate the integral 



corresponding to this scalar product. Repeating the proof of (|106|) given in || one 
can obtain a corresponding result in the case of the modified Bessel functions. The 
formula one gets turns out to be exactly of the form of (|106|) with A and B denoting 
now / or K. Thus we can study the scalar product similarly as we did in the case 
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s < 0. Let us evaluate ( |115[ ) with the aid of ( |106| ), the 'upper limit terms' give 
asymptotically 



V2 



V 2 ) 



c + (w)*c+(^)(e~™ - e w ) + c+(n)*c_(i;)(e-™ - e"™) 
+ c_(w)* C4 »(e™ - e 7 ™) + c_(M)*c_(w)(e" u - e" 7 ™) 



(116) 



Using the asymptotics ( ffiS|) we can see that the 'lower limit terms' behave as (a/2) l ( ±u± ^, 
just like in the case s < 0. Two of them vanishes if a — > 0. After a trick similar to 
( p.l4j ) the nonzero contribution of the two other terms is a sum of a sin A(u—v) / (u—v) 
and a cosA(w — v)/(u — v) term. In the limit A — > oo the first of them leads to 
a Dirac delta. The second one is not, consequently its coefficient must be equal to 
zero. This gives the condition 



c + {u)* c + {y) = c_(w)* c-(v). 



;ii7) 



On the other hand, (|116|) is a smooth function for u ^ v. Orthogonality requires 
that (|116|) must be equal to zero for any These two conditions together yield 



(uf + c-(u) m ] {c+(w)e™ + c-»e-™ - c + (v)e™ - c.(v)e~ wv } = 0. (118) 



[c+( 



( |118| ) holds for any u,v, u ^ v only if c + (u) + c_(it) = 0, which we wanted to prove 
(cf. @). 

We note that by using (|106|) and the asymptotics (|96| ) — (|1 02|) I v (z) and K v (z) 
prove to be not square integrable for any complex value of v. Another remark is 
that the eigenf unctions in (|63|) and flBlf) are real (cf. (pif)). 



B Completeness 

Here we prove the completeness of the orthogonal systems (|63|) and fl6"5|). Formulat- 
ing the completeness of a system (p3|) in the variable x reads 



with 



and 



S(xi,x 2 ) + I(xi,x 2 ) = S(xi - x 2 ) 

S(x 1 ,x 2 )= ^qJq{z\)J q {z 2 ) 

q£p+2Z+ 



I(x 1 ,x 2 ) 



du u 



o 2v / 2sinh 



7TU 



e- ie *Wj iu {zi)+e ie » {u) J- iu {zi 
e- ie * {u) J lu (z 2 )+e w ^J_ lu (z 2 ) 



x 



where 



Z\ = ke ^ , z 2 = ke V2. 



(119) 
(120) 



121^ 



(122) 
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The convergence of the infinite sum ( |120| ) is guaranteed by the asymptotics ( |103| ) . 
First let us prove (|119|) in the case when p = 2. We consider the integral (|121|) 



between and A, A = N + 1/2, iV G Z + . The substitution v := iu transforms the 
integral to 

[Mzi) + J-u(zi)] [Ju{z 2 ) + J- V {z 2 )\ . (123) 



k i2^2 sm7rz/ 

As the integrand of ( |123j ) is invariant under the transformation v — > —v ( |123|) can 



be written as 

1 r iA dv v 



2 J-iA i2 v / 2si 



smirv 



[Ju{zi) + J-u{zi)} [J v {z 2 ) + J- V {z 2 )\ . (124) 



Now we change the contour of this integral to a half circle (denoted by C) starting 
from the point — iA, running in the half plane Rez/ > of the complex z/-plane 
and ending at iA. The difference of (|124|) and this new integral can be expressed 
by the residues of the poles of the integrand lying in the region bordered by the 
two contours. As the J-s behave analytically, poles arise only from shi7rz/, at the 
values v — 1, 2, 3, . . . N (at v — vj sin(7rz/) is not singular). The contribution of 
the residues is 

2ni . , 

E Res « (125) 



AV2 



n=l 



where, by using (0), Res n turns out to be equal to J n {zi)J n (z2) if n is even 
and is zero if n is odd. We can see that ( |125| ) is just the opposite of (|120|) in the 
limit N — > oo so from now on we have to prove that 

[Jv{zi) + J-v(zi)} [J V (Z 2 ) + J- U {Z 2 )\ = 



c iA^/2 smnv 

dv v r dv v 



L 



. t /- . JM)Uz2) + / —jf^ J- v { Zl )J- v {z*) + (126) 

c i Ay 2 sum v Jc i4y/ 2sm7rv 

[Ju(zi)J- v (z 2 ) + J- v {z 1 )J v {z 2 )\ 



Ic i4\/2 sin irv 

tends to 5(x\ — x 2 ) when A — > oo. 

We will perform the proof in three steps. We start by showing that the first and 
the second terms of the r.h.s. of ( |126|) are equal. The second step proves that these 
terms are zero in the limit A — >• oo. Thirdly it is shown that the third term tends to 
5(xi — x 2 ) as A — » oo. 

To see that the first and the second terms are equal let us make the substitution 
v — ► — v in the second term. The integrand of the resulting integral is the same as in 
the first term (up to a factor of —1), but the contour C is a half circle starting from 
iA, running through the half plane Rez/ < and arriving at — iA. Let us change 
the direction of C , this causes another factor of —1 in the integral. Now we can 
change this contour C" to C , thus we arrive at the first term, plus the contribution 
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of the residues coming from the poles lying between C" and C. Fortunately this 
contribution 

N 

E ReS n (127) 
n=-N 

is zero because, as a consequence of fl9lD , Res n = — Res_ n . 

In the second step we are interested in the A — > oo behaviour of the integral 

(U ' U -J v {z,)J v {z 2 ). (128) 



ic smnu 

Writing v in the form of Aexp(iip) this integral can be expressed as 

I iA 2 e 2iip 

dip ^ Jxeiv{zi)J he i9{z 2 ). (129) 

s sm(7rAe* ¥ ') 

It will be enough if we show that the integral of the integrand's modulus tends to 
zero. 

We will use asymptotic expressions for analysing the A — > oo behaviour. By 
doing so some care will be needed. That's why we divide the domain of integration 
[-f, f] into three parts: A : = [-f,-f], B := [-f,|] and C := [f,|]. It can be 
proved easily that for a (p ^ and a large enough A 

|sin[7rAe i¥, ]| « i e - Asin M. (130) 

Let us apply ( |130|) in the domains A and C . Combining it with (|103|) the modulus 
of the integrand of ( |129| ) is asymptotically equal to 

^ e -2A[(lnA-C)cos V +(f-M)sinM] 5 ( 131 ) 
7T 

where C, — \ ln(ziz 2 /4) + 1 is a quantity independent of ip and A. We can see that 
the integral of (|131|) on A or C gives the same result. Thus we will consider this 
integral only on C, for example. 

If tp G [vr/4, 7r/2] then the inequalities cost/? > 1 — (2/ir)(p and (7r/2— yj) sin > 
hold, helping us to give an upper estimate of ( |131| ) 



^ e -2A(inA-0(i-| (132) 

7T 



if In A > The integral of ( |132j ) on C can be calculated easily. The result is less 
then [4(lnA — C)] -1 > which is a quantity tending to zero if A — ► oo. 

In domain B we cannot use ( |130|) but here it is enough to work with the inequality 

|sin[7rAe iv ]| > 1. (133) 

For proving ( |133| ) it is not hard to show that | sin[7rAexp(zy9)]| takes its minimum 
in (p = as (p varies in B while A is fixed. (Remember that A = N + 1/2.) For the 
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Bessel functions ( |103| ) is applicable in B, too. Using ( |133| ) the asymptotics of the 
absolute value of the integrand in ( |129|) is not greater than 



_ — 2A[(lnA— C)cos^— ipainip] 

2n 



(134) 



If ip £ B then cos<^ > l/y/2 and ipsinip < ix/Ay/2 so an upper estimate of ( |134j) is 



A 
2^ 



-^A(lnA-C-f) 



(135) 



if In A > Q. Integrating ( J135 ) on B means simply a factor of n/2. We see that the 
value of the integral tends to zero in the limit A — ► oo. 

In the last step let us turn to the third term of the r.h.s. of ( |126| ). For great A-s 
this integral is asymptotically equal to 



dv v 



c i4\/2sin 



r (l + u r 



+ 



r(i - v) r(i + v) 



(136) 



(A is not an integer so ( |96"D is applicable). Knowing that V(l+v)V(l—v) = iru/ sin ttv 
and writing v as Aexp(i<p) we get the integral 



A 



47rV / 2 



dip e 



up 



6 V2 



{x 1 -x 2 )e i v 



+ e 



(137) 



(cf. ( |122| )). We determine this integral by expanding the exponentials in power 
series and integrating the terms individually. The sum of the results is 



A 00 

TlV 2 k=0 



+ q {x 1 -x 2/ 



2k 



{2k + l)\ 
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which is the power series of sin 775(^1 — ^2) / tt(^i — ^2)- This function tends to 
5(xi — x 2 ) if ^ — > 00, which we wanted to prove. 

The method of the proof of the case p = 2 can be applied in a straightforward way 
for any other values of p as well. If p = 1 then a factor of —1 appears at the first two 
terms of the r.h.s. of ( |126|) , which is of no significant importance in the proof. For the 
other possible values of p the remarkable difference is that poles come not only from 
sin(7rz/) but also from the analytic continuation of exp[+i8 p (u)}. Nevertheless, the 
sum of all residues appearing in the proof gives exactly — S(xi, x 2 ), which eliminates 
the first term in ( |119| ) — just as it happened in the case p = 2. Besides these extra 
poles the proof needs no serious modification. 

The proof given above also works for the orthogonal system ([J5|) of the case 
s > 0. Using (|93|) the residues turn out to be zero. What is left is just the same as 



we had in the case s < 0, p = 1, because of the common asymptotic behaviour of 
J v (z) and l v (z) (see (|103|) ). This way the case s > can be treated with the same 
tools as the s < one. 
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